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We compute the power spectra in the cosmic microwave background and cold dark matter (CDM) 
fluctuations seeded by strings, using the largest string simulations performed so far to evaluate 
the two-point functions of their stress energy tensor. We find that local strings differ from global 
defects in that the scalar components of the stress-energy tensor dominate over vector and tensor 
components. This result has far reaching consequences. We And that cosmic strings exhibit a single 
Doppler peak of acceptable height at high t. They also seem to have a less severe bias problem than 
global defects, although the CDM power spectrum in the “standard” cosmology (flat geometry, zero 
cosmological constant, 5% baryonic component) is the wrong shape to fit large scale structure data. 


Recent times have witnessed unprecedented progress in 
mapping the CMB temperature anisotropy and the large 
scale structure (LSS) of the Universe. The prospect of 
fast improving data has forced theorists to new standards 
of precision in computing observable quantities. The new 
standards have been met in theories based on cosmic in¬ 
flation 1^. Topological defect scenarios have been 
more challenging. In these theories, as the Universe cools 
down, high temperature symmetries are spontaneously 
broken. Remnants of the unbroken phase, called topo¬ 
logical defects, may survive the transition, and later seed 
fluctuations in the CMB and LSS. The defect evolution is 
highly non linear, thereby complicating the computation 
of these fluctuations. 

Last year saw a number of computational break¬ 
throughs in defect theories, partly related to improve¬ 
ments in computer technology. Most strikingly the 
method described in showed how one could glean 
from defect simulations all the information required to 
compute accurately CMB and LSS power spectra. This 
method was applied to theories based on global symme¬ 
tries. Work on cosmic strings associated with gauged (Ur 
local) symmetries appeared at about the same time [|[|| , 
but making use of rather different methods. 

In this letter we report on a calculation of the local cos¬ 
mic string power spectrum, using the method of ap¬ 
plied directly to local string simulations. In this method 
the simulations are used uniquely for evaluating the two 
point functions (known as unequal time correlators, or 
UETCs) of the defects’ stress-energy tensor. UETCs are 
all that is required for computing CMB and LSS power 
spectra. Furthermore, they are constrained by require¬ 
ments of self-similarity (or scaling) and causality, which 
enable us to radically extend the dynamical range of sim¬ 
ulations, a fact central to the success of the method. 

We believe that our work has significant advantages 
over |^|j. In Q string simulations are used directly as 
sources for the cosmological perturbations. As the au¬ 
thors point out, this means that one is severely limited in 
dynamic range by the string simulation itself. The UETC 
method allows us to cover the full dynamic range required 
for CMB and LSS computations. In Q one made use of 


an analytical model for strings, first proposed by one of 
the authors in . Although the model has been shown 
to approximate some of the UETCs quite well 0i , our 
direct use of string simulations is clearly an improvement. 
We show elsewhere how the model misses some key 
features found in simulations. 

Local strings have an extra complication over global 
defects, which stems from the fact that we are unable to 
simulate the underlying field theory. Instead, we approx¬ 
imate the true dynamics with line-like relativistic strings. 
This is thought to be reasonable for the large scale prop¬ 
erties of the stress-energy tensor, but we do not have a 
good understanding of how the string network loses en¬ 
ergy in order to maintain scaling. In any case, one must 
conserve the total energy momentum tensor, and so one 
is forced to make assumptions about which cosmological 
fluids pick up this deficit. It is often assumed that all the 
strings’ energy and momentum is radiated into gravita¬ 
tional waves, approximated by a relativistic fluid. This is 
by no means certain, and it may well be that the energy 
and momentum is transferred to particles , and hence 
to the baryon, photon and CDM components. 

We explore these possibilities in this paper, and one 
of the main results is that the matter power spectrum 
is very sensitive to the assumptions made about string 
decay. In particular, it is possible to reduce the bias at 
100 Mpc scale to 1.6, which runs against the current 
orthodoxy, that defects necessarily have a large bias at 
this scale. The shape of the CDM power spectrum is 
still glaringly different from the data [^. Regardless 
of assumptions made on string decay products, we see 
a fairly distinctive peak in the CMB power spectrum, 
differing from |^,|| and from global defect theories, albeit 
with no secondary oscillations as expected for most active 
sources. 

We proceed to describe in detail our calculation. The 
unequal time correlators are defined as 

(0^,.(k, r)0*^(k, t')) = r, r') (1) 

where 0^^ is the stress energy tensor, k is the wavevec- 
tor, and r and r' are any two (conformal) times. The 
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UETCs determine all other 2 point functions, most no¬ 
tably CMB and LSS power spectra CV and P{k). Realis¬ 
tic UETCs have to be measured from defect simulations, 
although analytical modelling [^t|J^] gives insights into 
the observed forms. 

We performed flat space cosmic string simulations us¬ 
ing the algorithm described in [Q. We used a previously 
developed code ^ implementing this algorithm. Flat 
space codes achieve great efficiency and accuracy by ne¬ 
glecting the effect of Hubble damping on the network, 
and by restricting the strings to lie on a cubic lattice. 
They are thought to give a good quantitative picture of 
a real string network on large scales, although the overall 
string density is probably too high. Another weakness is 
that we cannot model the reduction in the string density 
at the radiation-matter transition. However, this is only 
a 20% effect in the comoving energy density j^ . 

We performed simulations in 128^, 256^, 450^ boxes, 
with a cut-off on the loop size of four links. To evalu¬ 
ate the UETCs from the simulations we selected times 
in the range 0.17V < t < N/A, where N is the box size, 
when we were sure that the string network was scaling, 
and when boundary effects are still excluded by causal¬ 
ity. For each of these times we compute and Fast Fourier 
Transform the string stress-energy tensor Qfj,v We then 
decompose the 0^jy(k, r) modes into scalar, vector, and 
tensor (SVT) components (eg. |^). Isotropy guarantees 
that the only non-vanishing correlators involve compo¬ 
nents with the same transformation properties, and as 
one can show there are only 14 independent UETCs. 
We compute them by cross-correlating a target time in 
the middle of our time range with all other times, and 
averaging over several runs. 

Scaling and causalilty impose powerful constraints on 
the correlators [^, and can be used to extend the dy¬ 
namical range of the simulations. Scaling implies that 
C^iv,ai 3 {k,T,T') = c^y^ai 3 {kT,kT')/yfrr'. The scaling func¬ 
tions c^v^ai 3 (kT,kT') can be found from simulations, al¬ 
though they are noisy in the fcr, kr' « 0 region. Causal¬ 
ity may be used to reduce the noise, as it implies that 
the real space correlators of the fluctuating part of 
must be strictly zero for r > t + t' . Real space correla¬ 
tors C{r, T, t') may then be truncated at r > r -|- t', and 
inverted into C{k,T,T'). This procedure completes our 
knowledge of the functions at kr ft! 0. 

We will report details in |j^, but a striking feature of 
our results is the dominance of ©oo over all other compo¬ 
nents. The string anisotropic stresses are in the predicted 
[p^ ratios |0'®p : : |0^P of 3 : 2 : 4, as fcr —> 0. 

However |0ooP ^ |0'^P) and so scalars dominate over 
vectors and tensors. Also the energy density power spec¬ 
trum rises from a white noise tail at A:r ft 0 into a peak 
at fcr ft 20, after which it falls off. Subhorizon modes are 
therefore of great importance. 

These features consistently appeared for all box sizes, 
and are independent of the cutoff size imposed on the 
loops. In 1^ we show that scalar dominance arises be¬ 
cause the other correlators are suppressed by powers of 


the mean square velocity of the string network, equal 
to about 0.36 (l^, and by geometric factors. 

The UETCs c^y^a^ikT^kr') may be diagonalised ||] 
and written as 

Cf,y,ai3{kT, kr') = Y^>^‘^"^vl^},{kT)v^^ji{kT') (2) 

i 

where are eigenvalues. In general, defects are inco¬ 
herent sources for perturbations [p^ , which means that 
this matrix does not factorize into the product of two vec¬ 
tors v^y{kT)vafi{kT'). Standard codes solving for CMB 
and LSS power spectra assume coherence. However we 
see that an incoherent source may be represented as 
an incoherent sum of coherent sources. We may there¬ 
fore feed each eigenmode into standard codes to find the 
cf'’ and (k) associated with each mode. The series 

and Y provide convergent approx¬ 

imations to the power spectra. 

The response of radiation, neutrinos, CDM, and 
baryons to coherent sources may be computed with a 
Boltzmann code (see |]I| for formalism). A popular and 
fast implementation is CMBFAST |j^, which is accurate 
to about 1%. Sudden recombination approximation 
codes are even faster, but only achieve about 5% accu¬ 
racy. We experimented with an implementation of the 
sudden recombination approximation, a full Boltzmann 
code, and CMBFAST. We found that for all active per¬ 
turbations tested our sudden recombination code never 
differs from the full Boltzmann code by more than 10% 
in Ct and 5% in P{k). We have also reproduced the re¬ 
sults in with this code. Since the uncertainties in the 
string UETCs lead to much larger errors, we felt that a 
sudden recombination approximation was good enough. 

A significant difference between current codes for sim¬ 
ulating local cosmic strings and global defects is the 
former’s lack of energy and momentum conservation. 
Long strings loose energy and momentum to small loops 
which are excised from the simulation. These can also 
be thought of as representing the decay of the long 
strings into gravitational radiation or high energy par¬ 
ticles. This feature introduces two novelties in the cal¬ 
culation. Firstly one cannot measure a reduced number 
of defect correlators (typically 3 scalars, 1 vector, 1 ten¬ 
sor land determine the others by energy conservation (as 
in [^). Instead we must compute all 14 correlators, and 
from them infer the long string violations of energy and 
momentum conservation. 

Secondly we must model the real physics of the de¬ 
cay products of the long strings. To this end one can 
introduce an extra fluid, specified by 2 scalar equations 
of state (eg. = w^©^, and H'^'^ = 0), and a vector 
equation of state (H'^^ = 0). For gravitational radia¬ 
tion = 1/3. For a fluid of loops ft u^/3, with v 
the rms centre of mass velocity. We explore the range 
0 < <1/3, although = 0 is probably unphysi¬ 

cal, as it is difficult to envisage that all the energy of the 
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strings could end up as the mass energy of non-relativistic 
particles. 

Another possibility is that strings decay into very high 
energy particles which must scatter and eventu¬ 

ally thermalise with the background fluids. This process 
entails transfer of energy and momentum to radiation, 
baryons, and CDM. In such scenarios, in addition to be¬ 
ing active perturbations, strings would also seed entropy 
fluctuations We shall explore all these possibilities. 

String decay products are clearly the most uncertain 
aspect of cosmic string theory. By measuring the full 14 
UETC associated with long strings, we assume nothing 
about decay products when extracting information from 
simulations (unlike where long strings and decay prod¬ 
ucts are modelled together). The simulations will then 
also place constraints upon the decay products. 

In Fig. H we plot \/\l[l -|- l)C'f/27r], setting the Hub¬ 
ble constant to iJo = 50 Km sec“^ Mpc“^, the baryon 
fraction to fib = 0.05, and assuming a flat geometry, 
no cosmological constant, 3 massless neutrinos, standard 
recombination, and cold dark matter. We superimpose 
also current experimental points. The most interesting 
feature is the presence of a reasonably high Doppler peak 
at £ = 400 — 600, following a pronouncedly tilted large 
angle plateau. This feature sets local strings apart from 
global defects. It puts them in a better shape to face the 
current data. 

The CMB power spectrum is relatively insensitive to 
the equation of state of the extra fluid. We have plotted 
results for = 1/3,0.1,0.01. Dumping some energy 
into CDM has negligible effect. Small dumps into baryon 
and radiation fluids, on the contrary, boost the Doppler 
peak very strongly. We plotted the effect of dumping 5% 
of the energy into the radiation fluid. 

The LSS power spectra on the other hand is strongly 
dependent on . In Fig. || we plotted the CDM power 
spectrum P(fc) together with experimental points as in 
(P^ . The normalization has been fixed by COBE data 
points. We see that the peak of the spectrum is always at 
smaller scales than standard CDM predictions, or obser¬ 
vations. However the overall normalization of the spec¬ 
trum increases considerably as vj^ decreases. 

The CDM rms fluctuation in 8 ft,“^Mpc spheres is erg = 
.42, .61,1.8 for = 1/3,0.1,0.01. Hence relativistic 
decay products match well the observed erg « 0.5. On 
the other hand in 100 /i“^Mpc spheres one requires bias 
5ioo = o’foo^/^ioo = 4.9, 3.7,1.6 to match observations. 

Energy dumps into radiation have no effect on the 
CDM power spectrum. However if there is energy trans¬ 
fer into CDM or baryons, even with vs^ = 1/3, the CDM 
power spectrum is highly enhanced. We plot the result of 
a 5% transfer into CDM and a 20% transfer into baryons 
(with = 1/3) for which feioo = 2.0,1.5. 

Hence in our calculations local strings have a bias prob¬ 
lem at 100 /i“^Mpc, although its magnitude is not as 
great as found in Q. It depends sensitively on the de¬ 
cay products, being reduced if the strings have a chan¬ 
nel into non-relativistic particles, or if there is some en- 



FIG. 1. The CMB power spectra predicted by cos¬ 
mic strings decaying into loop and radiation fluids with 
= 1/3, 0.1, 0.01, 0. We have plotted (£(£ -I- l)Cf/27r)^^^ in 
nK, and snperposed several experimental points. The higher 
curve corresponding to = 1/3 shows what happens if 5% 
of the energy goes into the radiation fluid. 

ergy transfer into the baryon and CDM fluid. The main 
problem with strings in an D = 1, Qb = 0.05, Ha = 0 
CDM Universe is that the shape of P{k) never seems to 
match observations. This may not be the case if one in¬ 
troduces Hot Dark Matter, thereby free streaming the 
smaller scales, or an open Universe, for which the scale 
of the horizon at equality is different. We reserve for 0] 
exploring the space of cosmological parameters. 

We performed a large number of checks on our results. 
We found fast convergence with the box size. The COBE 
normalized Ci vary by less than 2% as we go from 128^ to 
256^, and by less than a percent from 256^ to 450^. The 
normalization itself varies significantly. Gfj, decreases by 
25% as we go from 128^ to 256^, but hardly changes from 
256^ to 450^. For a 256^ box a COBE normalization at 
£ = 5 produces Gfi « 1.0 x 10“® (with = 1/3). Due 
to the large tilt this normalization changes considerably 
with the value of £ where the fit is made. The COBE 
normalized LSS power spectrum changes by less than a 
percent with the box size. 

Another useful test consists of changing the discretiza¬ 
tion (or weighting) applied to fcr, in the expansion §)■ 
If the eigenmode expansion is indeed working one should 
observe changes in the spectra corresponding to individ¬ 
ual eigenmodes, but not in the totals. We found that 
equal weighting in kr and in log(fcr) lead to final spectra 
in excellent agreement, with the latter displaying much 
faster convergence. 

The results we have obtained are consistent with pre¬ 
vious arguments on scalar, vector and tensor modes in 
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FIG. 2. The power spectrum in CDM fluctuations for cos¬ 
mic strings, with = 0.01,0.1,1/3. We plotted also the 
standard CDM scenario prediction and points inferred by Pea¬ 
cock and Dodds from galaxy surveys. The top 2 = 1/3 

curves correspond to a 5% transfer into CDM, and a 20% 
transfer into baryons (top). 

defect theories. In 0. rigorous arguments on the ra¬ 
tios |0'®p : for modes at fcr « 0 were 

derived. It was then showed how these translated into 
ratios Cf : Cj : Cj, under certain conditions. Two 
of the conditions were the subdominance of modes in¬ 
side the horizon [kr > 5), and that anisotropic stresses 
should have a similar amplitude to the energy density. As 
we have pointed out, local strings violate both these con¬ 
ditions. Hence, although we have observed the predicted 
ratios for the anisotropic stresses, the argument need not 
apply to Cl- We also checked our CMB code by using the 
UETCs of [|| as sources, and were able to reproduced the 
results. The conclusion is that the CMB and LSS pre¬ 
dictions for local strings and global defects are different 
because their UETCs are indeed qualitatively different. 

If we take « 1/3 our results are close to those of 
(a bias at 100 /i“^Mpc of 4.9 instead of 5.4; a higher 
Doppler peak). In Q the defect energy-momentum ten¬ 
sor is modelled as a gas of randomly oriented straight 
string segments, with random velocities, whose length 
and number density depend on time in the correct way 
to obtain scaling. In we develop further this analyt¬ 
ical model and show how the original model may miss 
some key features found in simulations. Hence the dis¬ 
crepancy found is not altogether surprising. Also in ||] 
the extra energy-conserving fluid is relativistic and non¬ 
interacting. That some of our results are quite different 
is explained by our widening the range of possibilities for 
this fluid. 

We stress that the results in Q assume a rather differ¬ 


ent background cosmology (iJo = 80 Km sec“^ Mpc“^, 
and fib = 0.02). In we shall report results for these 
parameters. In the range f = 100 — 300 we observe a C/ 
shape similar to |^, but with a higher amplitude. 

In summary, we have computed the CMB and LSS 
power spectra for local cosmic strings, using extensive 
flat space string simulations to model the sources. We 
have explored the consequences of relaxing previous as¬ 
sumptions about the decay products of the strings. We 
find that the 100 h~^ Mpc bias problem and the absence 
of a Doppler peak, thought to be generic features of de¬ 
fects, may not be as severe for local strings as they are for 
global defects. It appears that CMB and LSS power spec¬ 
tra depend on the details of the defect considered, and 
more seriously in the case of local strings, on the physics 
of the transfer of energy and momentum to matter and 
radiation. In an Einstein-de Sitter CDM Universe, with 
fit, = 0.05 and Hq = 50 km s“^ Mpc“^, the shape of the 
CDM power spectrum cannot be made to fit the data 0] 
even with our relaxed assumptions. Other cosmological 
parameters remain to be explored . 
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